Robust nodal superconductivity induced by isovalent doping in 

Ba(Fei_^.Ru^)2As2 and BaFe2(Asi_^.P^.)2 
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We present the ultra-low-temperature heat transport study of iron-based superconductors 
Ba(Fei_a;Rua;)2As2 and BaFe2(Asi_j;P2,)2. For optimally doped Ba(Feo.64R.uo.36)2As2, a large resid- 
ual linear term kq/T at zero field and a yfli dependence of ko{H)/T are observed, which provide 
strong evidences for nodes in the superconducting gap. This result demonstrates that the isovalent 
Ru doping can also induce nodal superconductivity, as P does in BaFe2(Aso.67Po.33)2- Furthermore, 
in underdoped Ba(Feo.77Ruo.23)2As2 and heavily underdoped BaFe2(Aso.82Po.i8)2, kq/T manifests 
similar nodal behavior, which shows the robustness of nodal superconductivity in the underdoped 
regime and puts constraint on theoretical models. 

PACS numbers: 74.70.Xa, 74.25.fc, 74.20.Rp 



Since the discovery of high- Tc superconductivity in 
iron pnictides [l|, the electronic pairing mechanism 
has been a central issue 0. One key to understand it 
is to clarify the symmetry and structure of the super- 
conducting gap [3]. However, even for the most studied 
(Ba,Sr,Ca,Eu)Fe2As2 (122) system, the situation is still 
fairly complex 

Near optimal doping, for both hole- and electron-doped 
122 compounds, the angle-resolved photon emission spec- 
troscopy (ARPES) experiments clearly demonstrated 
multiple nodeless superconducting gaps Q , which was 
further supported by bulk measurements such as ther- 
mal conductivity [7H9| . On the overdoped side, nodal su- 
perconductivit y w as found in the extremely hole-doped 
KFe2As2 [13, 

while strongly anisotropic gap 9], 
or isotropic gaps with significantly different magnitudes 
[l^ [l3| were suggested in the heavily electron-doped 
Ba(Fei_xCoa;)2As2. On the underdoped side, recent heat 
transport measurements claimed possible nodes in the su- 
perconducting gap of hole-doped Bai_a;Ka;Fe2As2 with 
X < 0.16 14|, in contrast to the nodeless gaps found in 
electron-doped Ba(Fei_3;Co^)2As2 
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Intriguingly, nodal superconductivity was also found 
in optimally doped BaFe2(Aso.67Po.33)2 (?c = 30 K) 
[lil Il6| . in which the superconductivity is induced by 
isovalent P doping. The ARPES experiments ha ve g iven 
conflicting results on the position of the nodes |l7l. fisf. 
Moreover, previously LaFePO (Tc ^ 6 K) displays clear 
nodal behavior 1^21|, and recently there is penetration 
depth evidence for nodes in the superconducting gap of 
LiFeP {Tc - 4.5 K) The nodal superconductiv- 

ity in these P-doped compounds are very striking, which 
raises the puzzling question why the P doping is so special 
in iron-based superconductors. The theoretical explana- 
tions of this puzzle are far from consensus (23l - [26| . 

A recent proposal is that the nodal state in iron- 



FIG. 1: (Color online). Isovalent doping in the Fe2As2 slabs 
of BaFe2As2, by substituting As with P, or Fe with Ru. Both 
ways can induce superconductivity, and result in similar phase 
diagrams. 



pnictide superconductors, except for KFe2As2, is induced 
when the pnictogen height hpn from the iron plane de- 
creases below a threshold value of '-^ 1.33 A :22]. Ac- 
cording to this proposal, there may exist a transition 
from nodal to nodeless state tuned by hpn, for exam- 
ple, in underdoped BaFe2(Asi_a;Pa;)2. Therefore, it is 
important to investigate the doping evolution of the su- 
perconducting gap structure in BaFe2(Asi_a;Pa;)2. In an- 
other aspect, since isovalent substituting Fe with Ru in 
BaFe2As2, as shown in Fig. 1, can also decrease hpn and 
result in similar phase diagram 27 - |29| , it will be very in- 
teresting to check whether the gap of Ba(Fei_a;Ru2:)2As2 
has nodes. 

In this Letter, we report the demon- 
stration of nodal superconductivity in opti- 
mally doped Ba(Feo.64Ruo.36)2As2, underdoped 
Ba(Feo.77Ruo.23)2As2, and heavily underdoped 
BaFe2(Aso.82Po.i8)2 by thermal conductivity mea- 
surements down to 50 mK. Our finding of nodal gap in 
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FIG. 2: (Color online). In-plane resistivity of 

Ba(Feo.64Ruo.36)2As2, Ba(Feo.77Ruo.23)2As2, and 

BaFe2(Aso.82Po.i8)2 single crystals. The low-temperature 
superconducting transitions are shown in the inset. Defined 
by p = 0, the transition temperatures Tc = 20, 17, and 4 
K are obtained, therefore these three samples are named as 
OP20K, UD17K, and UD4K, respectively. The solid line is a 
fit of the data between 30 and 90 K to p{T) = po + AT" for 
the OP20K sample. 



Ba(Feo.64Ruo.36)2As2 suggests a common origin of the 
nodal superconductivity induced by isovalent P and Ru 
doping. The nodal gap in Ba(Feo.77Ruo.23)2As2 and 
BaFe2(Aso.82Po.i8)2 shows no transition from nodal to 
nodeless state in the underdoped regime. 

Single crystals of Ba(Fei_a;Rua;)2As2 and 
BaFe2(Asi_:EPa;)2 were grown according to the methods 
described in Refs. [sO, |3l|. The Ru and P doping 
levels were determined by energy dispersive X-ray 
spectroscopy. The sample was cleaved to a rectangular 
shape with typical dimensions 1.50 x 0.7 mm^ in the 
a6-plane, and 40 to 80 ^m in c-axis. In-plane thermal 
conductivity was measured in a dilution refrigerator, 
using a standard four-wire steady-state method with 
two Ru02 chip thermometers, calibrated in situ against 
a reference Ru02 thermometer. Magnetic fields were 
applied along the c-axis. To ensure a homogeneous field 
distribution in the samples, all fields were applied at 
temperature above Tc. 

Fig. 2 shows the in-plane resistivity p{T) 
of Ba(Feo.64Ruo.36)2As2, Ba(Feo.77Ruo.23)2As2, and 
BaFe2(Aso.82Po.i8)2 single crystals. The transition tem- 
peratures defined by p = are Tc = 20, 17, and 4 
K, therefore we name these three samples as OP20K, 
UD17K, and UD4K, respectively. One can see clear 
resistivity anomalies in UD4K and UD17K, but not 
in OP20K, which manifest the gradual suppression of 
spin-density- wave (SDW) order upon P or Ru doping 
[ii, m, [ag. For OP20K, the resistivity data between 
30 and 90 K are fitted to p{T) = po + AT"", which gives a 




FIG. 3: (Color online). Low-temperature in-plane thermal 
conductivity of Ba(Feo.64Ruo.36)2 As2, Ba(Feo.77Ruo.23)2As2, 
and BaFe2(Aso.82Po.i8)2 in zero and magnetic fields applied 
along the c-axis. The solid lines are k/T = a -I- feT fit to all 
the curves, respectively. The dash lines are the normal-state 
Wiedemann- Franz law expectation Lo/po, with Lq the Lorenz 
number 2.45 x 10"* WfiK"^ and normal-state po = 44, 115, 
and 74 pficm, respectively. 



residual resistivity po = 43.7 ± 0.1 pQcm and n = 1.31 ± 
0.1. Such a non- Fermi-liquid temperature dependence of 
p{T) is similar to that observed in BaFe2(Asi_2;Pa;)2 near 
optimal doping, which may reflect the presence of anti- 
ferromagnetic spin fluctuations near a quantum critical 



point [32 1 



The resistivity of these samples were also measured in 
magnetic fields, the highest up to 14.5 T, in order to 
determine their upper critical field Hc2 and normal-state 
Po- For OP20K, UD17K, and UD4K, we estimate i/c2 = 
23, 19, and 5 T, and normal-state po = 44, 115, and 74 
pficm, respectively. 

Fig. 3 shows the temperature dependence of the in- 
plane thermal conductivity for OP20K, UD17K, and 
UD4K in zero and magnetic fields, plotted as k/T vs 
T. All the curves are roughly linear, as previously ob- 
served in BaFe1.9Nio.1As2 t8|], KFe2As2 llC|, and over- 
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doped Ba(Fei_a;Co2:)2As2 single crystals 0, There- 
fore we fit all the curves to k/T = a + bT°'~^ with a fixed 
to 2. The two terms aT and bT" represent contributions 
from electrons and phonons, respectively. Here we only 
focus on the electronic term. 

For OP20K in zero field, the fitting gives a residual 
linear term kq/T = 0.266 ± 0.002 mW K^^ cm"!. This 
value is more than 40% of the normal-state Wiedemann- 
Franz law expectation knq/T — Lq/pq — 0.56 mW 
cm~^, with Lq the Lorenz number 2.45 x 10~® 
W17K~^ and normal-state po = 44 /zficm. For optimally 
doped BaFe2(Aso.67Po.33)2 single crystal, similar value 
of kq/T « 0.25 mW cm"-'^ was obtained, which is 
about 30% of its normal-state kno/T [l^. The signif- 
icant Ko/r of Ba(Feo.64R.uo.36)2As2 in zero field is at- 
tributed to nodal quasiparticles, which is a strong evi- 
dence for nodes in the superconducting gap [ssl ]. 

With decreasing doping level, for UD17K and UD4K, 
Kq/T = 35 ± 1 and 69 ± 1 K'^ cm"! are obtained, 
as seen in Figs. 3(b) and 3(c). These values are about 
17% and 22% of their normal-state knq/T, respectively. 
We note Hashimoto et al. already mentioned that the 
superconducting gap in overdo ped SrFe2(Aso.6Po.4)2 and 
BaFe2(Aso.36Po.64)2 has nodes [2^. Therefore, our obser- 
vation of significant kq/T in the underdoped regime, par- 
ticularly in the heavily underdoped BaFe2(Aso.82Po.i8)2, 
further shows the robustness of nodal superconductivity 
against doping in P- and Ru-doped 122 iron pnictides. 

The field dependence of kq /T can provide further sup- 
port for the gap nodes In Fig. 4, the normalized 
{ko/T)/{kno/T) of OP20K, UD17K, and UD4K are plot- 
ted as a function of H/Hc2- For UD4K, k/T saturates 
above = 3 T, as seen in Fig. 3(c), which is determined 
as its bulk Hc2- For OP20K and UD17K, we use their 
Hc2 estimated from resistivity measurements. To choose 
a slightly different bulk Hc2 does not affect our discussion 
on the field dependence of kq/T. Similar data of an over- 
doped d-wave cuprate superconductor Tl-2201 [3J|, and 
BaFe2(Aso.67Po.33)2 are also plotted for comparison. 

For a nodal superconductor such as Tl-2201 in mag- 
netic field, delocalized states exist out the vortex cores 
and dominate the heat transport in the vortex state, 
in contrast to the s-wave superconductor. At low field, 
the Doppler shift due to superfluid flow around the vor- 
tices will yield an H^^^ growth in quasiparticle density of 
states (the Volovik effect ^35^), thus the H^^'^ field depen- 
dence of Kq/T. From Fig. 4, the behavior of kq{H)/T in 
OP20K, UD17K, and UD4K clearly mimics that in Tl- 
2201 and BaFe2(Aso.67Po.33)2- In the inset of Fig. 4, the 
KoiH)/T of OP20K, UD17K, and UD4K obey the H^/^ 
dependence at low field, which supports the existence of 
gap nodes. 

To our knowledge, so far there are five iron-based 
superconductors displaying nodal superconductivity, in- 
cluding KFe2As2 [lO|, lUj, underdoped Bai_a:Ka;Fe2As2 
(x < 0.16) 0, BaFe2(Asi_,P,)2 LaFePO [H- 



1.2 



1.0 



-i- 0.4 



0.0 



TI2201 

Ba(Fe Ru ) As 

^ 0.64 0.36'2 2 

Ba(Fe Ru ) As 

^ 0.77 0.23'2 2 



BaFe (As P ) 

2^ 0.82 0.18'2 




0.2 0.4 0.6 0.8 



0.0 



0.2 



0.4 0.6 

H/H ^ 



0.8 



1.0 



FIG. 4: (Color online). Normalized residual linear term 
Ko/T of Ba(Feo.64Ruo.36)2As2, Ba(Feo.77Ruo.23)2As2, and 
BaFe2(Aso.82Po.i8)2 as a function of H/Hc2- Similar data 
of an overdoped d-wave cuprate superconductor Tl-2201 [s^ ]. 
and BaFe2(Aso.67Po.33)2 [la] are also shown for comparison. 
The behaviors of k.q{H)/T in OP20K, UD17K, and UD4K 
clearly mimic that in Tl-2201 and BaFe2(Aso.67Po.33)2. In- 
set: the same data of OP20K, UD17K, UD4K, and Tl-2201 
plotted against {H / Hci)^^^ ■ The lines represent the H^^^ de- 
pendence. 



2l[, and LiFeP Here we only consider the "in-plane 
nodes", not counting the "c-axis nodes" in underdoped 
and overdoped Ba(Fei_a;Coa:)2As2 as suggested by c-axis 
heat transport experiments [36| . For the extremely hole- 
doped KFe2As2, the nodal gap may have d-wave symme- 
try, and result from the direct intra-pocket interaction via 
antiferroma gnet ic fiuctuations, due to the lack of electron 
pockets [l^, 37, [sl]. For underdoped Bai_a;Ka;Fe2As2, 
it is still not clear how the superconducting gap trans- 
forms from nodeless to nodal at x « 0.16 [ij]. The 
rest three compounds, BaFe2(Asi_a;Pa;)2, LaFePO, and 
LiFeP, have stimulated various interpretations of the ef- 
fect of isovalent P doping [23l - E6l |. and may represent a 
peculiar superconducting mechanism. 

Our new finding of nodal superconductivity in 
Ba(Fei_a;Rua;)2As2 reveals the similarity between the iso- 
valently Ru- and P-doped iron pnictides. In this sense, 
the P doping is not that special for inducing nodal su- 
perconductivity now, and the puzzle of P doping in iron- 
based superconductors has been partially unwrapped. 
What next one needs to do is to find out the common ori- 
gin for the nodal superconductivity in these isovalently 
doped iron pnictides. 

Due to the smaller size of P ion than As ion, one com- 
mon structural feature of the P-doped compounds is the 
decrease of pnictogen height hpn and increase of As-Fe- 
As angle [H, IS Ii3| • The substitution of larger Ru ion 
for Fe ion in Ba(Fei_a;Rua;)2As2 results in the increase 
of a lattice parameter and decrease of c lattice param- 
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eter, thus the decrease of pnictogen height and increase 
of As-Fe-As angle [28|. Therefore, both the P and Ru 
dopants cause the same trend of structure change in iron 
arsenides. 

With such structure change, the Fermi surface (FS) 
evolution upon isovlant P and Ru doping is rather del- 
icate. The main feature, hole pockets around Bril- 
louin zone (BZ) center and electron pockets around 
BZ corners, remains in LaFePO |41|, B aFe2(Asi-a;Pa;)2 



ITllialail, Ba(Fei_^Ru:,)2As2 |42|,|43|, and LiFeP |44|. 
For LaFePO, Kuroki et al. have attributed the low- Tc 
nodal pairing to the lack of Fermi surface 7 around (tt, tt) 
in the unfolded Brillouin zone, due to the low pnictogen 
height [11]. For BaFe2(Asi_a;Pa;)2, Suzuki et al. have 
proposed three-dimensional nodal structure in the largely 
warped hole Fermi surface and no nodes on the electron 
Fermi surface [2^. This is supported by recent ARPES 
experiments, which found nodal gap in the expanded a 
hole pocket at A:^ = tt in BaFe2(Aso.7Po.3)2 PLS|, how- 
ever, it conflicts with earlier ARPES results which have 
constrained the nodes on the electron pockets . Since 
ARPES experiments did not find significant changes in 
the shape of the Fermi surface or in the Fermi velocity 
over a wide range of doping levels in Ba(Fei_a;Rua;)2As2, 
Dhaka et al. speculated that its superconducting mecha- 
nism relies on magnetic dilution which leads to the reduc- 
tion of the effective Stoner enhancement [i^. In LiFeP, 
the middle hole pocket has significantly lower mass en- 
hancement than the other pockets, which implies that 
the electron-hole scatter rate is suppress for this pocket 
and may result in the lower Tc and nodal gap [i^ . 

While the clues for nodal superconductivity are not 
very clear from the FS topology except for KFe2As2, 
Hashimoto et al. gathered the available data for the low- 
energy quasiparticle excitations in several iron-pnictide 
superconductors, and suggested that there is a thresh- 
old value of hpn 1.33 A, below which all the su- 



perconductors exhibit nodal superconducting state [22 1. 
If this is the case, there may exist a transition from 
nodal to nodeless state tuned by /ip„ in underdoped 
Ba(Fei_a;Rua;)2As2 and BaFe2(Asi_a;P2;)2. To test this 
idea, we estimate hp,, = 1.317, 1.333, 1.340 A for OP20K, 
UD17K, UD4K from the roughly linear increase of /ip„ 
with decreasing Ru or P doping 2^ 3^ . 

One can see that both hpn of UD17K and UD4K are 
slightly larger than the proposed threshold value 1.33 
A. In particular, the /ip„ of UD4K is comparable to 
that of overdoped Ba.(Fen.s9Con.i 1 )2 As2, which is a node- 
less superconductor [22]. Since our thermal conductivity 
data suggest UD17K and UD4K are nodal superconduc- 
tors, hpn should not be considered as the only parame- 
ter for tuning between nodeless and nodal superconduct- 
ing states. By saying this, we do not deny its impor- 
tance, since hpn of the underdoped Ba(Fei_a;Rua;)2As2 
and BaFe2(Asi_a;Pa;)2 are still very close to the threshold 
value 1.33 A. More careful considerations of the struc- 



tural parameters, FS topology, and local interactions are 
needed to clarify the origin of the nodal superconductiv- 
ity in isovalently doped iron pnictides. 

In summary, we have measured the thermal conductiv- 
ity of Ba(Feo.64Ruo.36)2As2, Ba(Feo.77Ruo.23)2As2, and 
BaFe2(Aso.82Po.i8)2 single crystals down to 50 mK. A sig- 
nificant Ko/T at zero field and an H^/"^ field dependence 
of Ko{H)/T at low field give strong evidences for nodal 
superconductivity in all three compounds. Comparing 
with previous P-doped iron pnictides, our new finding 
suggest that the nodal superconductivity induced by iso- 
valent Ru and P doping may have the same origin. With 
decreasing doping level, nodal superconducting state per- 
sists robustly in heavily underdoped BaFe2(Aso.82Po.i8)2, 
suggesting that the hpn is not the only tuning parame- 
ter, thus putting constraint on theoretical models. Find- 
ing out the origin of these nodal superconducting states 
will be crucial for getting a complete electronic pairing 
mechanism in the iron-based high- Tc superconductors. 
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